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We present a numerical study of the self-affine profiles obtained from configurations of the g-state 
Potts (with q = 2, 3 and 7) and p = 10 clock models as well as from the occupation states for 
site-percolation on the square lattice. The first and second order static phase transitions of the 
Potts model are located by a sharp change in the value of the roughness exponent a characterizing 
those profiles. The low temperature phase of the Potts model corresponds to flat (a ~ 1) profiles, 
whereas its high temperature phase is associated to rough (a ~ 0.5) ones. For the p = 10 clock 
model, in addition to the flat (ferromagnetic) and rough (paramagnetic) profiles, an intermediate 
rough (0.5 < a < 1) phase - associated to a soft spin-wave one - is observed. Our results for 
the transition temperatures in the Potts and clock models are in agreement with the static values, 
showing that this approach is able to detect the phase transitions in these models directly from 
the spin configurations, without any reference to thermodynamical potentials, order parameters or 
response functions. Finally, we show that the roughness exponent a is insensitive to geometric 
critical phenomena. 
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I. INTRODUCTION 

In the past decade, the formation of rough surfaces un- 
der far-from-equilibrium conditions was a central theme 
in statistical physics. The application of self-affine frac- 
tals and scaling methods was essential to the progress 
that has been made towards the understanding of these 
non-equilibrium phenomena. Within this context, the 
standard tools used to describe various self-affine struc- 
tures observed in disordered surface growth are the 
roughness a and the growth /? exponents (l]-[| . The cen- 
tral goal of this approach is provide information about 
the correlations between fluctuations of a space and/or 
time varying property. Theoretical modeling of self-affine 
growth processes frequently used some of the models in- 
vestigated in critical phenomena, e.g., directed percola- 
tion, random field Ising Q and sine-Gordon models [Q. 
These efforts faced an old problem: thirty years after 
the renormalization group theory of critical phenomena, 
the quantitatively accurate prediction of the location and 
characteristics of phase transitions still constitute a chal- 
lenging and controversial question 0. This is partic- 



ularly true when one considers systems with random 
quenched disorder (e.g., random field Ising and Potts 
models, spin glasses and structural glasses ||), which 
exhibit non-trivial features such as extremely slow dy- 
namics, aging, ergodicity breaking, complex energy land- 
scapes, etc. 

On the other hand, the inverse problem, i.e., using the 
roughness exponents to study the main features of the 
phase diagram of equilibrium spin models, has not been 
much explored up to now. In 1997, de Sales et al. Q| 
mapped cellular automata (CA) configurations on solid- 
on-solid like profiles and used the roughness exponent a 
to classify the elementary Wolfram CA rules. Later, they 
also shown that this exponent can be used to detect the 
frozen-active transition in the one dimensional Domany- 
Kinzel CA (DKCA) |8| without any reference to order 
parameters or response functions. As we mentioned be- 
fore, beyond the roughness exponent a, the growth expo- 
nent (3 is another critical index used to describe rough- 
ening processes in the surface growth context. Recently, 
Atman and Moreira || determined the exponent (3 for 
the growth process generated by the spatiotemporal pat- 
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terns of the DKCA. The value of f3 exhibits a cusp at 
the frozen-active frontier and, if one observes the differ- 
ence configuration between two DKCA replicas, also at 
the active-chaotic critical frontier. The advantage of this 
method to find the phase diagram of the DKCA is that it 
is not necessary to wait for the system "thermalizes" , a 
process which often dispends a lot of computational time. 

In this paper we extended the roughness exponent 
analysis to other standard models of statistical mechan- 
ics. Specifically, we study the g-state Potts model (with 
q = 2,3 and 7), the simplest locally interacting statisti- 
cal model exhibiting both first and second order static 
phase transitions. We study also the p = 10 clock model 
for which a Kosterlitz-Thouless type phase transition is 
observed. Finally, since any connectivity problem can 
be studied by starting with pure random percolation and 
then adding interactions, we apply the roughness method 
to random-site percolation. 

In section 2 we describe the models and define the map- 
ping between spin configurations and walk profiles. In 
section 3, using the mapping between spin states and pro- 
files, we characterize the phases of the random-site perco- 
lation problem, Potts and clock models by the roughness 
exponent a. Finally, we conclude and indicate future 
directions of this work in section 4. 



II. MODELS AND FORMALISM 

The g-state Potts ferromagnet consists of spin vari- 
ables (Tj which may take on q discrete values cr^ = 
0,1, .... (q — 1) and are coupled by the dimensionless 
Hamiltonian 



(1) 



where S(, ) is the Kronecker delta function. 

The p-state clock model is defined by the Hamiltonian 



f3H = tfVcOsi— (<Ti-(Tj)\ 



(2) 



in which each spin can assume p discrete values <x; = 

0.1 V I- 

In Eqs. ([!]) and (^) the sums in (ij) runs over the 
lattice sites and their nearest-neighbors, (3 — 1/ksT, T 
is the temperature, K = f3J and J > is the coupling 
constant. We simulated the Potts model with q = 2, 3 



and 7 and the clock model with p — 10 states, on square 
lattices of N = L x L sites imposing periodic boundary 
conditions. For updating the spins we use a sequential 
Monte Carlo heat-bath process. 

In random-site percolation, one randomly occupies a 
fraction p of the sites of a <i-dimensional lattice (<Ji = 0: 
empty site; a = 1: occupied site). When p is small, 
the pair connectedness length scale £ p is short, compa- 
rable to the lattice constant a. However, when p ap- 
proaches p c , there occur fluctuations in the character- 
istic size of clusters on all scales from a to £ p , which 
diverges as £ p ~ \p — p c \ ~ yp . Each feature of thermal 
critical phenomena has a corresponding analog in perco- 
lation, so that the percolation problem is called a geo- 
metric or connectivity critical phenomena. For site per- 
colation on the square lattice the critical probability is 
Pc = 0.59275 ±0.00003 |u|. 

In the present work we mainly focused on the numer- 
ical study of the self-affine profiles generated from the 
configurations in the ordered and disordered phases of 
the Potts and clock models and in the site-percolation 
problem. As shown in a previous work [g the spin 
states can be mapped on random walk-like profiles and 
the correlations present in them can be measured using 
the roughness exponent. The simplest method to gener- 
ate walk profiles from the spin configurations at a time 
t is a 1:1 mapping in which each spin state (Ji(t) is as- 
sociated to a step (to the right or to the left) of a one- 
dimensional walk. Specifically, to an unique spin con- 
figuration {(7i(i), 02(f), . . • , cx/v(i)} corresponds a spatial 
profile {hi(t), /12(f), . . . , ft,jv(f)}, given by the sequence of 
the walker displacements hi after i unit steps defined as 



(3) 



tjj — (q — l)/2 for the q = 3 and q = 7 Potts 
= cry - 5 (p-j = <7j - 4) if a-j < 4 (if cr, > 4) for 



where pj = 
models, pj 

the p = 10 clock model, and pj = 2a j — 1 for the Ising 
(q = 2) model and site-percolation. 

After obtained the profiles by this mappings, the 
roughness exponent a was calculated by determining the 
average standard deviation of parts of the profiles with 
various scales e. At site i, in the scale e, the rms displace- 
ment fluctuation is given by: 



Wi(N,e,t) 



\ 



1 



2e+ 1 



(4) 



j=i-e 
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with 

i+e 

The roughness in the scale e is given by 
i * 

l^e.t) = . (6) 

i=l 

The roughness W(e) can distinguish two possible types 
of profiles. If it is random or even exhibits a finite correla- 
tion length extending up to a characteristic range (such 
as in a Markov chain), then W ~ e 1 / 2 as in a normal 
random walk. In contrast, if the profile has infinitely 
long-range correlations (no characteristic length), then 
its roughness will be described by a power law scaling 
such as 

W(e) ~ e a , (7) 

with a 1/2. The case a > 1/2 implies that the profile 
presents persistent correlations, i.e., a given displacement 
sequence (increasing or decreasing) is likely to be close 
to another of the same type. On the other hand, profiles 
with a < 1/2 are anticorrelated, which means that dis- 
placement sequences containing a great fraction of steps 
to the right are more likely to alternate with another 
one in which steps to the left are predominant and vice 
versa. The exponent a is restricted to the interval [0, 1] 
and is related to the fractal dimension d of the profile by 
a = 2-d@,|. 

III. RESULTS 

Before discuss the results, we shall make a very brief 
review of the equilibrium phase diagrams of the g-state 
Potts and p = 10 clock models. The square lattice q- 
state Potts ferromagnet presents a second order phase 
transition for q < 5 and a first order one for q > 5 at 

the critical temperatures Tc^ = l/ln(l + Jq) (in units 
of J/ks)- The p-state clock model interpolates between 
the Ising (p = 2) and the XY (p — > oo) models. For p > 5 
one expects the emergence of a soft spin-wave phase be- 
tween the ordered, low temperature, and the disordered, 
high temperature phases For the p = 10 case, the 

spin-wave phase is limited by the transition temperatures 
T ( /= 10) ~ 0.24 and T^ =10) ~ 1.0 @ . 



In Fig. 1 we show typical walk profiles generated by 
spin configurations of the Ising model (q = 2). From 
data similar to those in Fig. 1 for different temperatures, 
we obtained the behavior of the roughness exponents a 
as shown in Fig. 2, for the Ising model with L = 64 
and 128. These results correspond to averages over typ- 
ically M — 100 random initial configurations taken af- 
ter thermalization. The roughness exponent exhibits an 
abrupt fall from a ~ 1.0 to a ~ 0.5 at the tempera- 
ture T q= 2 ~ 1.13 which agrees with the static critical 

value Tc q ~ 2 ^ = 1.134... As we increase the system size we 
can note that the change of a at T q= 2 becomes sharper, 
suggesting a step function for the a(T) curve at the ther- 
modynamical limit, and a high-temperature value of the 
exponent a tending to a = 0.5. Indeed, the inset in Fig. 
2 shows that the width AT = T\ — goes to zero as 
N increases, where AT is arbitrarily defined as the dif- 
ference between the observed upper temperature T\ for 
which a ~ 0.999 and the temperature Ti for which a has 
fallen to half its maximum (a ~ 0.75). We also observe 
that, in this same limit, the temperature T% tends to the 
exact critical temperature Ti 9_2) {T\{L = 64) = 1.075, 
Ti(L = 128) = 1.10 and T X (L = 256) = 1.125). 

In Fig. 3 we show similar behaviors for the q = 3 
and q = 7 Potts model with L = 128. The abrupt 
falls in the exponent a are located at the temperatures 
T g= 3 ~ 0.99 and T q= ^ ~ 0.77 which are also in good 
agreement with the static critical values T , j 9_3 - ) = 0.994... 
and Tc q ~ 7 ' 1 = 0.773... The values a ~ 1 observed in 
the low temperature phases correspond to flat profiles 
reflecting the existence of long-range order (magnetiza- 
tion). In contrast, the values a ~ 0.5 obtained in the 
high temperature phases characterize random walk pro- 
files, as expected for disordered spin configurations. Fi- 
nally, any qualitative difference is observed among the 
a(T) curves across the critical surfaces corresponding to 
second (q = 3) or first (q = 7) order phase transitions. 

In Fig. 4 we show a typical log-log plot of W x e for 
the Ising model, whose fitted slopes give the roughness 
exponent values. We observe, in general, the existence 
of two distinct linear portions in the curve, whose inter- 
section point defines a particular length scale e*. The 
value of e* seems to be a measure of the average size 
of the spin islands or magnetized micro-domains limited 
by the correlation length £(T). Indeed, at the ferromag- 
netic phase, for small scales (e < e*) there is resolution to 
see the local spin fluctuations around the smooth profile, 
which leads to 0.5 < a < 1.0. For large scales (e > e*) 
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the profile is flat and the values a ~ 1 reflect the long- 
range order. At high temperatures e* rapidly decreases 
to the lattice constant a = 1, which agrees with the fact 
that such crossover is not observed in a random pro- 
file. An additional support to the relationship between 
e* and £ is provided by Fig. 5 which shows that e* (for 
the fsing model) has a peak near the temperature T q —2. 
The height and sharpness of this peak increase with the 
system size. A finite size analysis for the Ising model 
shows that the maximum value of the crossover length 
e* (which occurs near the critical point) scales linearly 
with the system size N = L 2 , similarly to the correlation 
length £. In fact, Fig. 5 suggests that e* diverges at T c as 
e* ~ (T — T c )~^ for the Ising model. Moreover, the be- 
havior of e* near Tc q ~ 2 ^ is similar to the 2D Ising model 
correlation length, whose exact expression £ = l/(/3* — /3) 
is known for /3 < (3 C jl3| (see Fig. 5). Here, the dual in- 
verse temperature j3* is given by (e 13 — 1) X (e 13 — 1) = 2. 

A recent work by Kantelhardt et al. Jl4[] provides the 
strongest support to the close relation between the length 
scales e* and £ that we suggest here based on our limited 
numerical evidence. These authors demonstrated that 
the detrended fluctuation analysis (DFA), which we used 
here in zero-order, can detect crossovers in the observed 
long-range correlation behavior of data series. They an- 
alyzed artificial data with a crossover from long-range 
correlations (a > 0.5) for s < s x to uncorrelated behav- 
ior for s > s x or vice-versa. Their DFA results clearly 
revealed the crossover and provided estimated crossover 

(n) 

lengths s x always larger than the real s x by a system- 
atic deviation which increases with the detrending order 
n. Also, the estimated s x n ^ were less accurate for a close 
to 0.5. The support for our suggestion comes from the 
analogy between the artificial series generated by Fourier 
transforms studied in jlj] and our data series build from 
spin configurations in which only thermal correlations, 
extending up to the scale £ (the correlation length), are 
present. 

The existence of two distinct linear portions in the 
plots of W x e is also observed in the Potts model. In Fig. 
6 it is shown the behavior of e* as function of temperature 
for the Potts model with q = 3 and 7. We can note the 
sharp peaks around the temperatures X'i 9-3 '' and Ti 9-7 ' . 

In Fig. 7 we show typical walk profiles generated by 
spin configurations of the p — 10 clock model at three 
temperatures: T — 0.15, located in the ordered phase, 
T = 0.40, in the spin-wave phase and T = 1.60 in the 
disordered one. The roughness exponent a as a function 



of the temperature for the p = 10 clock model is shown in 
Fig. 8. It suggests the existence of three distinct phases, 
namely, a low temperature flat (a ~ 1) phase with 
long-range order and extending up to Ta p =w — 0.25, 
a high temperature, disordered and rough (a ~ 0.5) 
phase for T > Tb p =io — 1-0, and an intermediate, 
rough (0.80 < a < 0.90) phase on the temperature range 
Ta p =iq < T < Tbj,=iq. These transition temperatures 
are in good agreement with those obtained for the static 
transitions between the ferromagnetic, paramagnetic and 
soft spin- wave phases of this model |T3] . The behavior of 
e* as a function of temperature is shown in Fig. 9. e*(T) 
exhibits a sharp peak around Tb p =io, the critical tem- 
perature separating the paramagnetic and the spin- wave 
phases, as well as a sudden jump near Ta p =w, where 
the transition between the ferromagnetic and spin-wave 
phase occurs. But, the main feature of Fig. 9 is that, 
in the spin-wave phase (Ta p =io < T < Tg P =io), e* 
has a high and almost constant value as expected for a 
Kosterlitz-Thouless "critical" phase characterized by an 
infinite correlation length £ at all those temperatures || . 

Finally, the roughness exponent method applied to the 
random site-percolation problem results in a constant 
value a = 1/2 for the roughness exponent over the en- 
tire range of the probability p. This result is expected 
since in site-percolation a fraction p of the lattice sites is 
randomly occupied, generating random profiles in both 
phases. In addition, it is not observed a length scale e* 
in the log-log plots of W x e, as occurred in the magnetic 
models. 



IV. CONCLUSIONS 

In this study we have shown that spin configurations 
for the Potts (with q = 2, 3 and 7 states) andp =10 clock 
models exhibit distinct self-affine characteristics, mea- 
sured by the roughness exponent a. The low temperature 
phases of the Potts model correspond to flat (a ~ 1) pro- 
files, whereas the high temperature phases are associated 
to rough (a ~ 0.5) ones. For the p = 10 clock model, in 
addition to the flat (ferromagnetic) and rough (param- 
agnetic) profiles, an intermediate rough (0.5 < a < 1) 
phase, associated to a soft spin-wave one, is found. The 
transition temperatures between the different roughness 
regimes are in good agreement with the static critical 
temperatures of these models. Our results show that the 
roughness exponent method is able to detect equilibrium 
phase transitions and provides accurate numerical deter- 



4 



mination of the critical surfaces without any reference 
to thermodynamical potentials, order parameters or re- 
sponse functions. In contrast, the roughness exponent 
method applied to the random site-percolation problem 
results in a constant value a = 1/2 over the entire range 
of the probability p, showing that the roughness exponent 
a is insensitive to detect geometric critical phenomena. 



The same analysis can be applied to damage-spreading 
phase transitions by focusing on the difference configu- 
rations between two replicas of the system, as done by 
Atman and Moreira ||. Moreover, using the growth ex- 
ponent /3, we can determine the phase diagram of spin 
models without any thermalization, leading to an impres- 
sive gain in simulations speed. 



We are extending our simulations in order to reliably 
estimate the value of the exponent /i and compare this 
value with the known critical exponents v for the Ising 
model. We conjecture that the exponent fi controlling 
the divergence of e* is the same as the correlation length 
exponent v for this model. 
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FIG. 1. Walk profiles obtained from equilibrium spin con- 
figurations of the Ising model (q = 2) at T = 1.50 > T c (<7=2) 
and T = 1.12 < T^ q=2) (inset). 



Another problem we can address using this rough- 
ness method is the location of phase transitions in disor- 
dered models such as random field and Ising spin glass 
models. For spin glasses, exchange (parallel tempering) 
techniques greatly improve traditional Monte Carlo al- 
gorithms, but the sizes and temperatures accessible to 
simulations are still insufficient to clearly solve several 
important questions In particular the value of the 
critical temperature for the d — 2 ± J Edwards- Anderson 
model is controversial. 



Apart from the application of the roughness method to 
spin models, a central issue is the nature of the relation- 
ship between the correlations measured on the spin-state 
profiles by the a or (3 exponents and the traditional two- 
point correlation functions C(r). We intend to examine 
this question more carefully in a future work. 
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FIG. 2. The roughness exponent a as a function of tem- 
perature for the Ising model with system sizes L = 64 and 
L = 128. The arrow indicates the static critical temperature 
j,jg-2) rpj^ mge £ gj lows the behavior of the width AT (see 
text) with the inverse of the system size N, indicating that 
a(T) is a step function at the N — > oo limit. 
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FIG. 4. A typical log-log plot of W(e) x e used to de- 
termine the roughness exponent a characterizing a profile 
obtained from an Ising spin equilibrium configuration for 
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FIG. 3. The roughness exponent a as a function of tem- 
perature for the q = 3 and q = 7 Potts model with L — 128. 
For q — 3 the model exhibits a second order phase transition, 
and for q = 7 a first order one. The arrows indicate the static 
critical temperatures T^ -3 ^ and ri 9-7 ' . 
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FIG. 5. The characteristic length e* as a function of tem- 
perature T for the Ising model with system size L = 128. e* 
exhibits a maximum around the critical temperature of the 
model. The solid curve represents the exact behavior for the 
correlation length £ in the disordered phase. The arrow indi- 
cates the static critical temperature T< 
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FIG. 6. The characteristic length e* as a function of tem- 
perature T for the q = 3 and g = 5 Potts models with system 
size L = 128. 
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FIG. 8. The roughness exponent a as a function of temper- 
ature for the p — 10 clock model with system size L — 128. 
This model exhibits a Kosterlitz-Thouless type phase transi- 
tion at T ~ 0.24. 
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FIG. 7. Walk profiles obtained from equilibrium spin con- 
figurations of the p — 10 clock model. T = 0.15 is located 
in the ordered phase, T = 0.40 is in the spin-wave phase and 
T = 1.60 is in the disordered one. 
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FIG. 9. The characteristic length e* as a function of tem- 
perature T for the p = 10 clock model with system size 
L — 128. The arrows indicate the static critical temperatures 
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